Back-Reaction and the Trans-Planckian Problem of Inflation Revisited 
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It has recently been suggested that Planck scale physics may effect the evolution of cosmological 
fluctuations in the early stages of cosmological inflation in a non-trivial way, leading to an excited 
state for modes whose wavelength is super-Planck but sub-Hubble. In this case, the issue of how 
this excited state back-reacts on the background space-time arises. In fact, it has been suggested 
that such back-reaction effects may lead to tight constraints on the magnitude of possible deviations 
from the usual predictions of inflation. In this note we discuss some subtle aspects of this back- 
reaction issue and point out that rather than preventing inflation, the back-reaction of ultraviolet 
fluctuations may simply lead to a renormalization of the cosmological constant driving inflation. 

PACS numbers: 98.80.Cq. 



I. INTRODUCTION 

The most important success of the inflationary Uni- 
verse scenario |l| is that it provides a causal mechanism 
for the origin of the observed density fluctuations and mi- 
crowave background anisotropies (see also Refs. 0,0]). 
Key to this success is the fact that the physical wave- 
length corresponding to a fixed comoving scale is expo- 
nentially stretched during the period of inflation. Thus, 
provided that the period of inflation lasts sufficiently 
long, fluctuations on scales of cosmological interest today 
originate on sub-Hubble scales during inflation. Since 
inflation red-shifts all initial classical fluctuations, it is 
reasonable to assume that matter starts out in a quan- 
tum vacuum state (in the frame set by the background 
cosmology, see e.g. Ref. p| for a discussion). Each fluc- 
tuation mode thus starts out in its vacuum state at the 
time that the initial conditions are set up (e.g. the be- 
ginning of the period of inflation), it undergoes quantum 
vacuum oscillations while the wavelength is smaller than 
the Hubble radius, but freezes out when the wavelength 
equals the Hubble radius (see e.g. Ref. p for comprehen- 
sive reviews of the theory of cosmological fluctuations). 
Subsequently, the quantum state of the fluctuations un- 
dergoes squeezing on super-Hubble scales, and re-enters 
the Hubble radius during the post-inflationary Friedman- 
Lemaitre-Robertson- Walker (FLRW) phase as a highly 
squeezed and effectively classical state (see e.g. Ref. |7| 
for a discussion of the classicalization of the state) . 

However, as first pointed out in Ref. p|, this success of 
inflationary cosmology leads to an important conceptual 
problem, the trans- Planckian problem. Since the period 
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of inflation in typical scalar-ficld-drivcn inflationary mod- 
els is very long (see e.g. Ref. 9] for a review), the scales of 
cosmological interest today are not only sub-Hubble, but 
in fact sub-Planck at the beginning of inflation. Thus, the 
formalism used to calculate the evolution of fluctuations 
is in fact not justified. It is possible that the unknown 
trans-Planckian physics will lead to an evolution of the 
fluctuations on sub-Planckian scales which from the point 
of view of free scalar field theory coupled to General Rela- 
tivity looks non-adiabatic. In Ref. 0| (see also Ref. [Tll|). 
toy models for such an evolution were constructed mak- 
ing use of modified dispersion relations which were as- 
sumed to describe the physics on sub-Planckian scales. 
Since the time interval spent in the trans-Planckian do- 
main may depend on the wavelength, such models may 
lead to changes in the spectral index of the fluctuations. 

Subsequently, other approaches to the trans-Planckian 
problem were suggested, e.g. analysis based on 
space-space non-commutativity [T^ . space-time non- 
commutativity minimal trans-Planckian assump- 

tions (starting each mode in some vacuum state at the 
time when its wavelength equals the Planck length) 0, 
fill . These analysis typically give that trans-Planckian 
corrections to the predictions for cosmological fluctua- 
tions are proportional to (Hi n {/m c ) n where H ln t is the 
Hubble parameter during inflation, m c a new scale at 
which non standard physical effects show up and n a 
number which depends on the initial state assumed at the 
time of "creation" . On the other hand, analysis based on 
modified dispersion relations give a correction propor- 
tional to the time spent by the physical modes in the 
region where adiabaticity is violated, see e.g. for a 
recent review. 

Tanaka [13 and Starobinsky (see also Ref. |Tgjl 
have, however, raised an important concern regarding the 
possible amplitude of trans-Planckian corrections (see 
also Ref. 20] for unrelated concerns): if trans-Planckian 
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physics leads to an excited state for fluctuation modes 
on sub-Hubble but super-Planck scale during the period 
of inflation, the back-reaction of these excitations on the 
background must be considered. In the case where the 
trans-Planckian effects are modeled by a modified dis- 
persion relation, a simple estimate of the energy density 
carried in these ultraviolet modes 

(P)uv = / d 3 fcph ys Wphys(fcphys)™fcphy S , (I) 

where fc p h ys is the physical wavenumber, Hfe . a is the oc- 
cupation number, and w p h ys is the frequency of the mode, 
leads to the conclusion that (p) uv will exceed the back- 
ground density unless n-k phyB is smaller than (H ln {/m rl ) 2 
[see also Eq. (Ql above] , thus constraining the possible ef- 
fects of trans-Planckian physics on the spectrum of fluc- 
tuations (in the above, we have assumed that the ultra- 
violet cutoff is the usual Planck scale m pl which is not 
mandatory at all). 

In this article, we point out some subtleties with the 
above back-reaction argument which may change the con- 
clusions dramatically. If we assume, as is conventionally 
done in analyzing quantum fields in curved space-time, 
that the ultraviolet cutoff scale is time-independent in 
terms of physical length, then in an exponentially ex- 
panding background geometry, the contribution of ultra- 
violet (i.e. sub-Hubble) modes to the energy density is 
constant in time, as, in fact, follows directly from the 
time translation invariance of the physics. Moreover, the 
corresponding equation of state, due to the fact that the 
dispersion relation is modified, can strongly differ from 
that of ultra-relativistic particles and, as we demonstrate 
below, tends to that of the vacuum. Hence, our main 
conclusion is that, instead of preventing inflation, the ul- 
traviolet modes may in fact simply renormalize the value 
of the cosmological constant driving inflation. 

This article is organized as follows. In the next sec- 
tion, Sec. II, we describe the arguments that have been 
put forward to claim that there is a back-reaction prob- 
lem and we criticize them. Then, in Sec. Ill, we present 
an explicit calculation of the equation of state of a scalar 
field with a modified dispersion relation. We show that 
the ultraviolet modes possess an equation of state which 
is almost that of a cosmological constant. Finally, in 
Sec. IV, we point out problems with our approach, indi- 
cate directions for further investigations and present our 
general conclusions. 

II. THE BACK-REACTION PROBLEM 

As explained in the introduction, it is possible that 
the trans-Planckian effects affect the standard inflation- 
ary predictions. In this paper, for the sake of illustration, 
we model physics at very short scales by a non-linear 
dispersion relation w p h ys (fcphys)- For wave- numbers such 
that fcphys "C m c where m c is a new scale at which non- 
standard physical effects show up, the dispersion relation 



is linear for obvious phenomenological reasons. On the 
contrary, for modes such that fc p h ys n% c > the shape of 
^phys (fcphys) is a priori unknown. It is has been shown 
that a non-adiabatic evolution of the mode function in 
the trans-Planckian region necessarily implies a modifica- 
tion of the inflationary predictions, in particular a modifi- 
cation of the power spectrum. To be more precise, in cos- 
mology, the dispersion relation becomes time-dependent 
and equal to u> = aaj p h ys (k/a) (u> and k denote comov- 
ing frequency and wavenumber, respectively). Then, the 
Wentzel-Kramers-Brillouin (WKB) approximation is sat- 
isfied provided that \Q/uj 2 \ -C 1, where the quantity Q is 
defined by Q = 3(cj') 2 /(4cj 2 ) — w" '/(2u) (a prime stands 
for the derivative with respect to conformal time), see 
also Ref. |2l|. If the previous condition is worked out, 
then one sees that the WKB approximation is violated if 
<^phys < -ffinf and that corrections to the standard result 
can occur in this case. 

This conclusion has been criticized in Refs. [ItI ITsf 
and many reasons why a modification of the inflationary 
power spectrum would be unlikely have been provided in 
these articles. In the following, we will examine each of 
them. 

In Ref. [lsf. it has been claimed that if w p h y s (fcphys) 
is such that the WKB approximation is violated for 
fcphys > Tn c then there are no preferred initial condi- 
tions. This is certainly correct for the class of dispersion 
relations considered in Ref. , as discussed in Ref. [23| , 
but not true in general. An explicit counter-example has 
been provided in Ref. [2^ and is studied in the present 
paper. The corresponding dispersion relation is sketched 
in Fig. 2] On this plot, one notices that in the region 
where fc p h ys > vn c there is an interval of finite range, 
namely fc p h ys € [Ai, A2], in which the WKB approxima- 
tion is violated (i.e. w p hy S < H- m f). Important for our 
analysis is also the fact that for /c p hy S — ► +00 the adia- 
batic approximation is restored, which follows since we 
have w p hys > Hmf. In this latter regime, the adiabatic 
vacuum is obviously the preferred initial state. 

The other arguments involve the calculation of the 
stress-energy tensor and are as follows. In the standard 
scenario, the initial conditions are fixed for all wavenum- 
bers at some initial time. If the number of e-foldings of in- 
flation is greater than about 70, the physical wavelength 
of modes which are currently probed in cosmic microwave 
experiments is smaller than the Planck length at the be- 
ginning of inflation. One usually assumes that the evolu- 
tion starts out from the adiabatic (Bunch-Davis) vacuum. 
If the evolution is non-adiabatic in the region fc p hys ^* > 
then the state in the region H m t < fc p h ys < rn c will 
differ from the usual adiabatic vacuum. Therefore, if 
one concentrates only on what happens in the region 
-ffinf < fcphys < Tti c , the trans-Planckian effects boil down 
to a modification of the initial conditions. This last argu- 
ment has been used in Refs. 0,0] as follows. Roughly 
speaking, a non- vacuum state means a non- vanishing en- 
ergy density and there is now the danger that this dom- 
inates over the energy density of the inflationary back- 
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ground which is mt .H?^ . According to Refs. 0,^1], this 
is actually what happens unless the level of excitation of 
the initial state compared to the adiabatic vacuum is very 
small, leading to unmeasurably small trans-Planckian ef- 
fects of the spectrum of fluctuations. Since observational 
evidence seems to indicate that inflation is the correct 
theory of the very early universe, Refs. [H[ll conclude 
that trans-Planckian effects of significant importance are 
in fact not possible. 

In order to understand the above argument in more 
details, let us be more accurate about what has actually 
been done in Refs. [TtI Hsj. It is well-known that cos- 
mological perturbations (density fluctuations and grav- 
itational waves) can, in some contexts, be viewed as a 



free scalar field (p(r),x) on a time-dependent background 
space-time. In the case of gravitational waves, the corre- 
spondence is exact, for scalar metric fluctuations (density 
perturbations), the correspondence is only exact if the 
equation of state of the background is time-independent. 
In the general case, the squeezing factor for the density 
fluctuations is given not by the FLRW scale factor a(t), 
but by a function z(t) which depends both on the back- 
ground geometry and the background matter - for details 
see e.g. Then, the corresponding energy density and 
pressure are given by the mean values of the stress-energy 
tensor (!),„). In an excited state characterized by the 
mode distribution function n = n(k), one has 
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where /ifc is the rescaled Fourier amplitude, i.e. /!& = 
a{rj)ipk{rj) and normalized such that /ik — l/\2k. In 
the above, 77 denotes conformal time, and a prime the 
derivative with respect to 77. 

In the above expressions, the terms proportional to 
the factor 1/2 are divergent in the ultra-violet regime, 
i.e. fc — > +00, and represent the quantum vacuum con- 
tribution. This means that, in order to give sense to the 
above expressions, the stress energy tensor should be first 
properly renormalized, for instance by adiabatic regular- 
ization |24j |. In this paper, we will simply subtract the 
contribution of the quantum vacuum energy since this 
is what has been done in Refs. 0, 0|. Another justi- 
fication is that we are in fact mainly interested in the 
terms proportional to n(fc) which describe the contribu- 
tions originatin g fr om the excited quanta. 

In Refs. 0,13; considerations have been restricted 
to physical modes such that H ln { < fc p h ys < Tn c . In this 
region, the dispersion relation is linear and, since the 
WKB approximation is satisfied, the mode function can 
be written as 

/i fc = ^Le- lfc " + & e +ifc " , (4) 

where |afe| 2 — \(3k\ 2 = 1- Inserting this mode function 
into the vacuum expressions of the energy density and 
pressure, one finds 

(p)w = A~i ^W, (5) 

27r 2 a 4 J aHiD{ fc 

1 1 r am c Ah 



We see that the coefficient \f3k\ 2 represents the number of 
particles, n(fc) = \(3k\ 2 - This describes the modification of 
the standard initial conditions due to the trans-Planckian 
effects (let us remind that the usual adiabatic initial con- 
ditions correspond to = 1 and (3^ = 0). Then, a back- 
of-the-envelope calculation shows that (p) uv — m 4 |/3fc| 2 
where we have used the fact that, in de Sitter space-time, 
the coefficient (3k is scale-independent (time translation 
invariance). If we require that the energy density of the 
test scalar field be smaller than the background density, 
then it follows that 

A similar expression has been obtained in Ref. [I^ > with 
to c = m P1 and, crucially, H m { replaced with Hq, the 
present value of the Hubble parameter. It is clear that 
the constraint on (3k is completely different (and much 
more difficult to satisfy) if one uses H ~ 10~ 61 m pl in 
the above equation rather than -ffi n f — 10~ 5 m pl (we also 
notice that m c needs not be the Planck mass) . The rea- 
son for this difference is again (see above the discussion of 
the preferred initial conditions) that, in Ref. it was 
assumed that violation of the WKB approximation in the 
trans-Planckian region necessarily implies that tu(k) — > 
as fc — > +00, as for the dispersion relation envisaged in 
Ref. |2^. In this case, the adiabatic condition is vio- 
lated today for a range of trans-Planckian modes, and 
one should indeed replace Hi n { by H as done in Ref. 0] . 
This results in a very strong constraint on (3k- This was in 
fact the essence of the criticism made in Ref. [23| against 
the dispersion relation considered in Ref. |22j. However, 
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again, the argument does not apply for dispersion rela- 
tions of the type shown in Fig. ^ an d, therefore, is not 
true in general. The reason can be very easily understood 
from Fig. ^ Since, after inflation, the Hubble parame- 
ter decreases and since, at some point, its value becomes 
smaller than the minimum of the dispersion relation, the 
adiabatic condition is restored at late times for all modes, 
and particle production stops. Therefore, in this case, 
the calculation should be done with H-,„f and not with 
H as done in Ref. |18|. In Refs. 0,123, it has been 
shown that there is a window for which the modifica- 
tion of the power spectrum is not completely negligible 
and for which there is no back-reaction problem, see for 
instance the discussion after Eq. (68) in Ref. 

The other arguments presented against possible trans- 
Planckian modifications of the inflationary power spec- 
trum involve the calculation of the equation of state 
(which is characterized by the parameter Lu st = p/p)- 
The main argument is that the "dangerous" created par- 
ticles behave as a radiation field (or as ultra-relativistic 
particles) and, therefore, that their energy density scales 
as a" 4 and not as the vacuum, see Ref. .17]. We now 
explain why this line of reasoning is problematic. 

First of all, using the stress-energy tensor of a test 
scalar field is questionable since one wants in fact to cal- 
culate the equation of state of the cosmological pertur- 
bations, not that of a test scalar field. It has been shown 
in Refs. (2^, [2(| that the equation of state of the effec- 
tive stress-energy tensor of the cosmological perturba- 
tions differs on super-Hubble scales (which is a region in 
which the adiabatic condition is not valid and thus has 
similarities to the trans-Planckian interval where the adi- 
abatic condition is violated) from the equation of state 
of a test field. It is cj s t = — 1 instead of ui s t = —1/3. The 
difference is clearly of utmost importance in the present 
context. In fact, what should be done is to calculate the 
stress-energy tensor of cosmological perturbations (which 
is second order in the perturbed metric) in the case where 
the dispersion relation is modified. To our knowledge, 
this calculation has never been performed, and is very 
complicated. Therefore, this is beyond the scope of the 
present article. 

Keeping the previous point in mind, let us come back 
to the calculation of the stress-energy tensor of a test 
scalar field. As already mentioned in the introduction 



[and demonstrated explicitly in Eq. 



fact constant and does not scale as a -4 despite the fact 
that p/p = 1/3. The non-conservation of the energy- 
momentum tensor can be understood as follows. Using 
the expression of the energy density, it is easy to establish 
that 
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from the time-dependent limits of integration in Eq. 0. 
The first one comes from the upper limit while the sec- 
ond one originates from the lower limit. The correspond- 
ing physical interpretation is clear: due to the expan- 
sion of the background, there is a flow of modes coming 
from the trans-Planckian region and entering the region 
H m f < kphys < m c and there is also a flow of modes 
leaving the region H ln { < /c p hys < m c while they are be- 
coming super-Hubble modes. Eq. (JSJ is similar to Eq. (4) 
of Ref. 18] . The only difference is the absence in Ref. 
of the second term on the right-hand-side (describing the 
outgoing flow of modes). This term is necessarily present 
because the integral in Eq. cannot be computed with 
a vanishing lower integral since, for fc p h ys < -ffinf, the 
mode function is no longer given by Eq. but rather 
by [ik ~ cl(ji). However, the second term in JHJ) is clearly 
very small in comparison with the first one and, there- 
fore, can be safely neglected. 

Finally, maybe the most important reason why calcu- 
lating the back-reaction can be more subtle than pre- 
viously thought is the following. The conclusion that 
p/p = 1/3 is in fact obtained from an inconsistent 
procedure since it does not take into account the fact 
that the dispersion relation is modified [let us remind 
that the previous considerations are based on Eqs. 
and that have been obtained under the assumption 
that Wphys = fcphys]- I n the following, we shall study the 
equation of state of the ultraviolet terms 10 and in 
a toy model for trans-Planckian physics in which we can 
describe the excitation of the mode functions in the far 
ultraviolet range from well-defined vacuum initial condi- 
tions in a mathematically consistent way. We shall show 
that the fact that the dispersion relation is modified can 
change the equation of state, a conclusion also reached 
in Ref. H3 in a slightly different context. 

To summarize, the calculation of the back-reaction 
must be performed with the trans-Planckian corrections 
taken into account (i.e. in the present context with a 
modified dispersion relation) . It is clearly inconsistent to 
calculate the modified power spectrum with the trans- 
Planckian corrections on one hand and, on the other 
hand, to evaluate the corresponding back-reaction with- 
out these corrections. This can change, in a crucial way, 
the calculation of the energy density and/or the equation 
of state. 



III. THE TOY MODEL 
A. Description of the Dispersion Relation 

We model the trans-Planckian effect by means of the 
following non-standard dispersion relation [23J 

W phys (^phys) = ^phys ~~ 2&nfc phys + 2&i 2 fc phys ■ (9) 



The two terms in the right-hand-side of the above ex- 
pression, responsible for the non-conservation, originate 



This dispersion relation is chosen such that the modes 
evolve adiabatically for extremely high wavenumbers, 
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but, given an appropriate choice of the constants &n and 
&12, there is an intermediate region of wavenumbers in 
which the mode evolution is not adiabatic. We will start 
the modes in their adiabatic vacuum in the extreme ul- 
traviolet and calculate how they are excited during the 
phase in which the evolution violates the adiabaticity 
condition. 

If we introduce the new dimensionless coefficients a 
and (3 such that a = 2bnm 2 , and (3 = 2b±2m 4 ,, where m c 
is a new free energy scale to be specified later on, then 
the dispersion relation can be re-written as 



m c J 



phys 



phys 



•j8 



^phys 



(io) 

It is represented in Fig. ^ This relation is in fact char- 
acterized by one parameter, the "shape parameter" T 
defined by T = 3/3/a 2 . The derivative of the dispersion 

relation vanishes at kifi/m c — ^/a/(3/3)y/l =F \/l — T 
which shows that T < 1 and the requirement that the 
dispersion relation stays positive implies that T > 3/4. 
To summarize, one has 



Let us notice that H- m f crosses the dispersion relation 
three times only if Q 3 + R 2 < 0. This implies that i?i n f 
should be chosen such that 



H m i n < H in f < H 



max ) 



(16) 



with -ff m i n and H max given by the following expressions 
(which, obviously, only depend on the shape of the dis- 
persion relation) 
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Obviously, the scales ki and k% only depend on the shape 
of the dispersion relation, i.e. only on the parameters a 
and p. It is more convenient to express everything in 
terms of a and T. This gives 



Then the three solutions can be found explicitly since 
they are in fact solutions of a third order polynomial 
equation (more precisely, u> 2 h — H 2 af is a sixth order 
polynomial equation that can be reduced to a third order 
equation in the variable k 2 hys ). The three solutions can 
be written as 



m 2 c 



- T 



(12) 
(13) 



In this paper, for simplicity we restrict ourselves to 
the case where the background space-time is de Sitter, 
characterized by the constant Hubble parameter H lTL i- 
Then, fco , Ai and A2 are the two scales for which w p hy S = 
-f/inf, see Fig.^ They depend on the parameters a and 
/3 but, clearly, also on H m {/m c . In fact, their explicit 



expressions can easily be derived. For this purpose, 
us define the coefficients Q and R by 
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(21) 



where 6 = cos 1 (R/y / — Q 3 ). On can check that, if 3/4 < 
T < 1, then k < Ai < A 2 as required. 
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Let us now consider a scalar field the dispersion rela- 
tion of which is given by Eq. I|1(J[I ■ Then, as demonstrated 
in Ref . |23l | , the vacuum expectation value of the energy 
density and pressure are given by 



(P) 
(P) 



1 



Aw 2 a 4 



r+00 

/ dkk 2 
la 


a 2 


a' 


r+co 

/ dkk 2 
la 


a 2 





- u) 2 (k) \nk\ 



x 2 2 dw 2 2 , 1- 



(22) 
(23) 



^phys/ m C 



FIG. 1: Sketch of the dispersion used for the study of the toy model. The five scales ko, ki, fe, Ai and A2 are defined in the 
text. ifi n f is the Hubble parameter during inflation while Ho is the Hubble parameter today. Clearly, this is not a scaled figure 
since Ho should be much smaller than is represented on the plot. 



This stress-energy tensor is conserved, i.e. (p)' + 3H{p + p) — 0, with H = a' /a as shown explicitly in Ref. [2e| . 
The modification of the energy density has exactly the expected form while the modification of the pressure is more 
complicated, involving the derivative of the dispersion relation. One can easily checked that, for the linear dispersion 
relation, the above formulas reduce to the standard ones. 



B. Near ultra-violet region 



Let us now try to evaluate these expressions explicitly. If we are in a region where the WKB approximation holds, 
then the mode function can be written as 
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with |cv(fc)| 2 — |/3(fc)| 2 = 1 from the Wronskian normalization condition. Inserting this expression into the formula 
giving the energy density, one obtains 
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where we have used la*! 



\Pk\ 2 - In the above expression, the quantity 7 is defined as follows 
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In a situation where WKB is a good approximation, we 
have 7/w ~ i and the previous expression reduces to 
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where the domain of integration K, corresponds to the re- 
gion where the WKB approximation is valid. Note that 
in order to remove the two oscillatory terms, no proce- 
dure of time averaging is needed in contrast with what 
was done in Ref. [l7| . 

We now demonstrate that the energy density of 
Eq. H27JI is constant in time. Replacing the time- 
dependent comoving frequency a; (A;) by its expression in 
terms of the physical frequency, namely u> = auj p hy S (k / a) , 
one gets 



(P) = 
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(28) 



where we have used that is scale-independent in the 
case of a de Sitter background. We have also specified the 
domain of integration JC. It is of course crucial that this 
domain be defined in terms of physical wavenumbers. As 
announced before, (p) is time- independent. Of course, 
this does not imply that the corresponding equation of 
state is necessarily —1 because, since we consider only a 
limited range of wavenumbers, the energy density is not 
conserved. Only the total energy density, integrated over 
wavenumbers from to +00, is conserved. 



Let us now evaluate the pressure. Repeating the same 
steps as before, a long but straightforward calculation 
gives 



(?) 



f 



47r 2 a 4 



2uj 



dkk 2 { — 

2 2 duj 2 

3 d*2 ~ 



2,2dw 2 2 1 I |2 

3 fc dP ~ W +W 



\Pk\ 



-2i f w(fe,T)dr , ^ffe 



LU 

2u 



2 2 du 2 

3 dfc 2 ~ 
'2 2 dtu 2 

3 dP 



+ (7*) 2 



3 2i f u(k,r)&T 



(29) 



This time, in order to remove the oscillatory terms, we 
can take the time average of the previous expression, as 
done for the energy density in Ref. 17]. In the follow- 
ing we denote the corresponding double average by the 
symbol ((■■•)). This yields 



1 1 



3 27r 2 a 4 



dkk 2 



+ \P k \ 2 \u{k) 



dlnw 2 
' d In k 2 



(30) 

The expression for the pressure can be evaluated explic- 
itly. We use the fact that the coefficient 0k is scale- 
independent, thanks to the time translation invariance. 
Then 



«P» 



l m 4 |/3 fc | 
3 2tt 2 



Ai 



V 



fcp 



(31) 



where the function Viz) is defined by the following ex- 
pression 



V(z) = 



T In 



16a 2 T 2 

27 

32a 2 T 2 



with lo(z) = z 2 — az 4 + a 2 Tz 6 /3. 



(9 - 8T) + 2aTz 2 + -a 2 T 2 z 4 
o 



-aTz 2 



(32) 



From the above expression, the energy density can be 
evaluated very simply if one notices that 



dfcfc 2 ^(fc) ^ hl ^ = k 3 uj(k)\ic-3 [ dkk 2 oj{k). (33) 
K din A; 2 J K 



Then the equation of state oj s % = p/p can be expressed 
as 



(A x /m c 



(fc /m c 



m c / (Ai/to c ) — T 5 (k /m c ) 



1 



(34) 



A typical example is represented in Fig. [21 The strik- 
ing feature of this plot is that, for H. af /m c -C 1, the 
equation of state goes to one. Let us emphasize that this 
regime corresponds to a physical requirement that should 
be met if one wants to be in a realistic situation. In order 
to understand better how this happens, we perform the 
following perturbative treatment. 

If the Hubble parameter is small in comparison with 
the new scale m c , then the minimum of the dispersion 
relation is close to zero which in turn means that T ~ 
3/4. Admittedly, this is a fine-tuning of the shape of 
the dispersion relation. The above considerations suggest 
that that Taylor expansion in T — 3/4 can be performed. 
Then, one finds [see Eqs. and 



R 1 , 27 (H inf 

— , = — 1 H a 

4 \™c 



■OIT--1. (35) 
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FIG. 2: Left panel: Equation of state versus Hi n f/m c for a = 0.01 and T = 0.7003 (dotted line) computed according to 
Eq. JHll' The dashed line is an approximation of the equation of state, valid for small values of H[ n t /m a , and derived in 
Eq. gTJ. Right panel: same as left panel but with a = 0.0001 and T = 0.700003. 



The next step is to expand the cos -1 function which ap- 
pears in the paragraph following (|21|l with the above ar- 
gument. At this point, since we expect H m { -C m c , one 
can treat a(Hi n {/m c ) 2 as a small parameter and Taylor 
expand our expressions, e.g. Eq. I|35[l . in this parameter. 
However, this approximation will break down when H ln { 
becomes large, especially far from H m i n - This is true even 
if T = 3/4 since then a(H max /m c ) 2 = 8/27 ~ 0.3 which 
represents an error of ~ 30%. In the vicinity of H m i n the 
approximation is of course much better. Then, working 
at zeroth order in T — 3/4, it is easy to show that the 



angle 9 introduced in the paragraph after (|21|l is 




+0 (r - | J . (36) 

If one inserts the above expression into the formulas giv- 
ing kg and Ai, one obtains 




The first expression is expected since it says that at leading order fco — -ffinf- This is because the usual transition 
between sub and super-Hubble modes occurs in the region where the dispersion relation is almost linear. In the same 
manner, one has 




Finally, putting everything together, at leading order in T — 3/4 and in a(ffi n f/m c ) 2 , we obtain a simple equation 
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in the regime of interest, namely 




(41) 



At some time ^(k), the comoving mode enters the far 
ultra-violet region, i.e. when the physical wavelength 
equal A2, and fik and its first derivative must be matched 
to the above solution, which gives 



It is represented by a dashed line in Fig. [3 We are now 
in a position where the behavior of the equation of state 
can be understood better. In the limit T — » 3/4 and 
H in { — ► -ff m in, the equation of state goes to —1. Let us 
notice that, since we have H^ in = 8(T — 3/4)/(3a) + 
• • • , the limit H — > -ffmin corresponds, in this case, to 
Hi n { — ► 0. We conclude that, in the physical regime 
of interest, namely H/m c <C 1 the equation of state is 
extremely close to that of the vacuum. Far from H m i n 
the approximation used above breaks down as is apparent 
from Fig. When H m i — > H max , it is clear that k — > Ai 
and the expression giving the equation of state becomes 
ambiguous. From the plots, we see that w s t goes in fact 
to zero. Therefore, even if the Hubble constant is not 
small in comparison with m c , the equation of state to a t 
remains negative. 



; J^W ^(fc, r )dr 0.(7?) 



x \ 1 - 7 (fc,7yi) 



?J2 W 



o(t) 



dr , (43) 



where the quantity 7 has been defined previously, see 
Eq. I|26(l . Then, one can safely neglect the second term 
which is the decaying mode and express everything in 
terms of physical quantities, using that the scale factor 
can be written as a(r)) = —l/(H ln {rj). This gives 



A 2 



ui[k,rj 2 (k)] fc2 hys (r?) 



(44) 



C. Far ultra-violet region 

Let us now consider the far ultra-violet region. We fix 
the initial conditions in the region where the WKB ap- 
proximation holds by selecting positive frequency modes, 
which corresponds to the choice of the Bunch-Davies adi- 
abatic vacuum state for the field, 



1 



VMM) 



exp 



i I w(fc,r)dr 



(42) 



But, one has uj[k,rj 2 (k)] = a(?7 2 V ph ys[fc/a(?72)] = 
a(??2)wphys(A 2 ) = a(r) 2 )H m {. Since a[r) 2 (k)] = fc/A 2 , one 
finally arrives at 



(45) 



2fcff inf fc2 hys ( f? )- 
Then, it is straightforward to calculate the energy density 
and the pressure. This gives 



<P> 



1 



16^ 2 H iDt 



1 



Kim 2 



48tt 2 Hu 





•Ar 






(-)'- 


( A 







a 
2 



A, 
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2 



A, 



m„ 



Ai_ 




3 



A, 



A, 



m„ 



Ai 

m„ 



(46) 
(47) 



In the standard case where a = f3 = 0, one recovers that p/p = —1/3 as required for super-horizon modes (see e.g. 
|26|L Then, straightforward algebraic manipulations show that 



(P) 



-<P) 
-<P> 



2 Ah 



3 i?inf 



2 A2to 2 



3 



inf 



A, 



Ai 



A, 



i^phys (A 2 ) - ^phys (Al)] = ~(p) 



Ai 



A 5 



(48) 
(49) 



Therefore, in the far ultra-violet region, the equation of 
state is nothing but the vacuum equation of state, i.e. —1. 

Can we understand this result better? For this purpose 
let us consider Eqs. I|22|l again. The terms (/ifc/a)' vanish 



because fik oc a in the region under consideration. Then, 
the link between the pressure and the energy density can 
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be re-written as 



(P) = -(P) + 



1 



3 47r 2 a 4 



(50) 



Now, if |/ife| 2 scales as \p,k\ 2 cx l/k 3 as indicated by 
Eq. (|45|l which is the consequence of having matched the 
mode function in the far ultra-violet region to the initial 
Bunch-Davis vacuum, then the above expression can be 
re- written as 



(p) + (p) <* 



dfc- 



di^ 2 
dfc 



(51) 



Clearly, if the interval K. is such that the frequency is 
the same at its boundaries then we obtain the equation 
of state of the vacuum. This conclusion does not de- 
pend on the detailed shape of the dispersion relation in 
this region. It is also obvious that the calculations done 
previously for a specific dispersion relation are fully com- 
patible with the above considerations, in particular the 
difference between (p) and -(p) in Eq. (|48|l is given by a 
term equal to 2/3 times the difference between the square 
of the effective frequency at the boundaries of the far 
ultra-violet region. 

To conclude this section, one can estimate the equation 
of state of the created particles coming from the whole 
ultra-violet region. In the near ultra-violet region, one 
has 



(P) 



near — UV 



-l + O 



(P) 



near— UV j 



while in the far ultra-violet region 



(p) 



far-UV 



= -(p) 



far-UV 



(52) 



(53) 



In the "ultra- far" region (fcphys > A2), i.e. in the region 
where the initial condition are fixed, we have by defi- 
nition the adiabatic vacuum and, therefore, no created 
particles. Hence, we do not need to take into account 
this region. As a consequence, the equation of state of 
the whole ultra-violet region can be written as 



^st-UV 



_ (P> 



near— UV 



+ (p> 



far-UV 



(p)near-UV + (p) far-UV 



l + O 



(54) 

where we have used Eqs. (15211 and i|53[l . Since the total 
energy density 



(p)uv = (p)near-UV + (p) far-UV ■ 



(55) 



is constant in time, we see that the energy density of the 
created particles almost behave as a positive cosmologi- 
cal constant. The slight non-conservation is due to the 
fact that we have not taken into account the infra-red 
region which also contributes (let us recall that the to- 
tal energy- momentum tensor is conserved exactly). The 
same calculation performed with a linear dispersion re- 
lation would have led to the result w st _uv = 1/3- We 



have thus demonstrated explicitly that taking into ac- 
count the trans-Planckian corrections is important when 
one evaluates the back-reaction. 

We end this section with the following remark. It is 
important to keep in mind that the previous calculation 
is not the calculation of the equation of state cosmologi- 
cal perturbations. Therefore, one cannot claim that the 
equation of state of the cosmological perturbations with 
a modified dispersion relation is the one of a cosmological 
constant. What has been calculated is just the equation 
of state of a test field. Nevertheless, we have shown with 
the help of the previous toy model that any attempt to 
evaluate the equation of state of the cosmological per- 
turbations in a regime where the dispersion relation is 
modified must take into account these trans-Planckian 
corrections. 



IV. DISCUSSION AND CONCLUSIONS 

Now that we have determined the energy density and 
pressure of the ultraviolet modes, we will study their 
back-reaction on the background space-time and mat- 
ter. We work in the context of the toy model studied 
in the previous section. We assume that at early times 
all modes with wavenumber larger than A2 start out in 
their adiabatic vacuum. Thus, after subtraction of the 
quantum vacuum terms, these modes do not contribute 
to the energy-momentum tensor, and thus there will be 
no remaining ultraviolet divergences. 

We now follow a mode with fixed comoving wave- 
length which starts out deep in the ultraviolet region 
(&phys > A2) in its adiabatic vacuum. The mode will then 
spend a finite time interval in the intermediate frequency 
range Ai < /c p hy S < A2 during which the adiabaticity 
condition for mode evolution is violated, the state gets 
squeezed, and, from the point of view of the vacuum state 
for fcphys < Ai, a non- vanishing occupation number nk 
is generated. This occupation number remains constant 
when ko < fc p hy S < Ai, since in this frequency range the 
adiabaticity condition is restored. 

Given this setup, we consider the energy density puy of 
ultraviolet modes (modes with wavelength smaller than 
the Hubble radius) which is the sum of l|27|l and 14611 . 
see Eq. (|55|l . By time-translation invariance of the back- 
ground, both terms are independent of time in a de Sitter 
background in which Hi„f is constant. Thus, 



dt 



0. 



(56) 



This seems to indicate that the ultraviolet energy den- 
sity evolves like a cosmological constant, and its back- 
reaction, rather than preventing inflation, will simply 
lead to a renormalization of the cosmological constant. 

Naive intuition, namely treating the equation of state 
as radiative, i.e. p uv = p vv /3, would have led to a 
problem with this conclusion, namely an extreme non- 
conservation of the energy-momentum tensor of the ul- 
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traviolet modes, a non-conservation on the energy density 
scale of Af (in this context the difference between Ai and 
A 2 is not relevant). However, our analysis of the previ- 
ous section has shown that in fact the energy- momentum 
tensor of the ultraviolet modes is that of a cosmological 
constant, up to correction terms which are suppressed by 
a factor of H ln f/m c compared to the dominant terms. 

Let us take a look at the equations of back-reaction. 
Following the method discussed in Refs. [2^, |3(j in the 
context of the problem of the back-reaction of infrared 
cosmological fluctuations. The back-reaction effect of 
interest here is the fact that linear cosmological fluctu- 
ations effect the background metric and matter if one 
works out the Einstein equations to quadratic order in 
the amplitude of the primordial perturbations. To be 
specific, in the following we shall consider a homoge- 
neous, isotropic and spatially flat background, and will 
take matter to be a scalar field p with a quadratic poten- 
tial given by the scalar field mass m. For the fluctuations 
we parameterize the equation of state as 



(-1 + aK 



(57) 



where (by the above discussion) a is a positive constant 
expected to be of the order Hi n { /m c . 



Our starting point consists of taking the FLRW equa- 
tions 

H 2 = |p, ^ = -3H(p + p), (58) 

where k = 8Tr/m 2 r The presence of the fluctuations 
produces a back-reaction effect on the background which 
is quadratic in the amplitude of the fluctuations [2^, |26| 
and leads to a correction of both metric and matter, i.e. to 
corrections SH of the Hubble expansion rate and dip of 
the scalar field 



H = H lnf + SH, ip = p in{ + S<p . 



(59) 



where the subscripts "inf" stand for the quantities evalu- 
ated in the unperturbed background. The back-reaction 
of the linear fluctuations on the background is described 
by contributions pbr = P vv and pbr = P vv (and we use 
the ultraviolet energy densities and pressures as the back- 
reaction quantities) to the energy density and pressure. 
If we insert the back-reaction ansatz (|59f) into the equa- 
tions H58fl (written in the slow-roll approximation) and 
linearize in SH and Sip, we obtain the following equa- 
tions 



2H ini SH = - (p uv + m 2 <p ini 5tp) 



( d<A, 



V dt 



Sip 



m 2 ip- u 



/ dStp 
\~dT 



-3H hl {(p l]V +p uv ) 



d/\ 



dt 



(60) 



r 



We have assumed that the background model has a mas- 
sive potential, i.e. V{<p) = m 2 ip 2 /2. The back-reaction 
of infrared modes (modes with wavelength larger than 
the Hubble radius) was analyzed in |2(j (see also 
Ref. [U) and was shown to correspond to a negative 
cosmological constant whose absolute value increases in 
time (see Ref. 32] for resulting speculations on how this 
effect might be used to address the cosmological constant 
problem). Here, we will study the back-reaction effects of 
the ultraviolet modes for which dpbr/dt = dp uv /dt = 
and for which the equation of state is given by (|57|l . In 
this case, the linearized Klein-Gordon equation can be 
simplified to yield 



dS(p 
~df 



, r dyw , 



Vinf 



3aH„ 



9 Puv 



The solution of Eqs. (jHU and JHOl is 



Sp = 
SH = 



-3 

K p v 



\/!(— 

V 6 V m 



6 H 



(1 



inf 



at , 

3aH mf t) 



(61) 

(62) 
(63) 



As expected Sp and SH vanish if p 
33] and later analyzed in detail in [3 



= 0. As stressed in 
it is important to 



express the result in terms of physical observables instead 
of in terms of the non-measurable background time coor- 
dinate t. The obvious clock in our simple system is the 
scalar field ip = ip in f + Sp itself. After a simple calcula- 
tion it follows that the back-reaction effect SH measured 
in terms of ip is given by 



H in{ + SH= x j- 



mp - 



6 H 



inf 



(64) 



Thus, the effect of the ultraviolet back-reaction terms 
in a de Sitter background corresponds to a positive re- 
normalization of the cosmological constant. This result 
does not depend on the value of a. 

So far, we have shown that a large ultraviolet back- 
reaction does not prevent inflation, in contrast to naive 
expectations. Instead, it leads to a re-normalization of 
the cosmological constant. On the other hand, as is evi- 
dent from Eq. I|63|l . the back-reaction terms can lead to 
a faster rolling of the scalar field i^mf- However, for the 
small values of a of the order of Hi n {/m c indicated by 
our analysis, the increase in the rolling speed does not 
prevent a phase of inflation of sufficient length (i.e. hav- 
ing I Sp> I < Iv'inf I for a time interval t ~ H^ { ) as long as 
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n k < — , (65) 
to p , 

which is a less stringent constraint than the one derived 
by Tanaka, and still allows observable effects on cosmic 
microwave fluctuations from trans-Planckian physics. 

We have presented an attempt to study the effects 
of back-reaction on the trans-Planckian problem of in- 
flationary cosmology. The question initially posed in 
Refs. 17] and [3 is whether the back-reaction of a state 
which consists of excited modes during the inflationary 
phase on scales smaller than the Hubble radius will pre- 
vent inflation. We have seen that the analysis is more 
subtle than it initially appears from the above works. We 
have shown that the back-reaction of an excited state 
preserving time-translation invariance does not prevent 
inflation, but simply leads to a renormalization of the 
Hubble constant. However, back-reaction will lead to a 
slightly faster rolling of the scalar field. As long as the 
occupation numbers are smaller than H{ n { /m rl , the 
rolling will be consistent with the standard inflationary 
paradigm. Such occupation numbers can lead to observ- 
able effects on the cosmic microwave background. 

To render the analysis well defined, we have consid- 
ered a dispersion relation for which the violation of adi- 
abaticity is concentrated in a finite range of physical 
wavenumbers Ai < fc p h ys < A2. We have assumed that 
all modes start off in their adiabatic vacuum state when 
&phys > A2. They are squeezed while Ai < k p i iys < K%, 
and then emerge as excited states when fc p hy S < Ai. Let 
us also notice that by making the wavelength interval 
Ai < fcphys < A2 small, the energy density in the far ul- 
traviolet modes can always be made small compared to 
the cutoff energy density Af (or Af). Thus, m our ap- 
proach the Planck energy density problem recently dis- 
cussed in Ref. [3(| in an approach to quantum field theory 
on a growing lattice, in which the number of fundamen- 
tal field modes is increasing in time, does not arise (the 
analysis of Ref. [3(| finds that the continual creation of 



modes at the cutoff scale yields an energy density which 
is of Planck scale). In our analysis, the Hubert space 
of modes is time-independent, but as time proceeds an 
increasing subset of this space gets populated. Our mech- 
anism for continual excitation of new modes appears to 
be more smooth than the lattice approach of Ref. [3(j . 

There are several important deficiencies in our anal- 
ysis. First, we have used an ad-hoc regularization and 
renormalization prescription which consists of imposing 
effectively an abrupt cutoff in momentum space, and of 
subtracting the ground state energy of each field Fourier 
mode. This procedure is not covariant. It would be of 
interest to study our problem using a mathematically 
more rigorous regularization prescription, such as adi- 
abatic regularization (see e.g. Ref. [37( and references 
therein). 

Another serious concern is that we have considered 
matter fluctuations without taking into account the in- 
duced metric fluctuations. It is well known that the in- 
clusion of metric fluctuations leads to dramatically differ- 
ent results for super-Hubble-scale perturbations [2{| |26| . 
Thus, one might expect that the gravitational fluctua- 
tions could play an important role in the far ultravio- 
let region where w p h ys < -Hinf- However, at the present 
time we are not able to study this issue because the non- 
standard dispersion relation has been set up for the mat- 
ter sector only. It is a challenge for future research to 
include the presence of non-standard dispersion relations 
consistently in both the gravitational and matter sector. 
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